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Abstract: The OPERA Collaboration reported evidence for muonic neutrinos traveling 
slightly faster than light in vacuum. While waiting further checks from the experimen- 
tal community, here we aim at exploring some theoretical consequences of the hypothesis 
that muonic neutrinos are superluminal, considering in particular the tachyonic and the 
Coleman-Glashow cases. We show that a tachyonic interpretation is not only hardly recon- 
ciled with OPERA data on energy dependence, but that it clashes with neutrino production 
from pion and with neutrino oscillations. A Coleman-Glashow superluminal neutrino beam 
would also have problems with pion decay kinematics for the OPERA setup; it could be 
easily reconciled with SN1987a data, but then it would be very problematic to account for 
neutrino oscillations. 

PACS numbers: 14.60.St 



I. INTRODUCTION 

Recently the OPERA Collaboration [l] reported an early arrival time for CNGS muon 
neutrinos with respect to the one expected assuming neutrinos to travel at the speed of 
light in vacuum c. The relative difference of the velocity of the muon neutrinos v with 
respect to light quoted by OPERA is: 

= (2.48 + 0.28{stat) + 0.30(sys)) x 10"^ , (1) 
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for a mean energy of the neutrino beam of 17 GeV. 

Similar hints, but with much less significance, were also reported for muon neutrino 
beams produced at Fermilab. Dealing with energies peaked at 3 GeV, the MINOS Collab- 
oration m found in 2007 that {v - c)/c = (5.1 + 3.9) x 10"^. In 1979, a bound on the relative 
velocity of the muon with respect to muon neutrinos (with energies from 30 to 200 GeV) 
was also extracted: \v - v^i\lv^i < 4 x 10"^ [i3j|. 

While urging the experimental community to further check and debate on these results, 
in particular the most recent ones, it is worth to explore which theoretical consequences 
would follow from the hypothesis that the muon neutrino is a superluminal particle. 

This clearly requires a deep modification of the Standard Model (SM) of particle 
physics, that assumes particles to be subluminal. The energy and momentum of a sublu- 

minal particle of mass m and velocity v are: E - mc^/ ^jl — ^, p = mv/ ^jl — ^. They are 
related by the dispersion relation = p^c^ + m^c^, where p = \p\, and the deviation from 
the speed of light is: 



c + v 



2\2 



(2) 



Here we focus on two possibilities to account for a superluminal particle: the tachyon 
IHl and the Coleman-Glashow particle Il6j[7|. 

The usual expressions for energy and momentum can be extended to the region v > c 
provided we substitute in the numerator m — > im, where m is a real number. For such 

a particle, the energy and momentum are thus E = fhc^/ - 1, p = mv/ — 1, and 
satisfy the dispersion relation 

£2 = p2c2 - m2c4 . (3) 
A tachyonic particle of mass mc^ and energy E then travels faster than c by an amount 



V -c c mc ' 



c + v 



(4) 



The deviation from the speed of light for a tachyon is thus opposite with respect to the 
one of a subluminal particle with the same mass, eq. Q. In both cases, for energies 
much larger than the mass, the particle speed v approaches c. We display in fig. [Tjthe 
relative speed of the tachyon with respect to light as a function of the tachyon mass and 
for selected energy values. 

The proposal that the neutrino could be a tachyon dates back to 1985 ||4|. In the light 
of the recent results, it is worth to consider this hypothesis as a possible explanation of 
the OPERA data. As we are going to discuss in section |ll| this interpretation is very 
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problematic for various reasons: not only the deviation from c would depend on energy, 
but it would not even be possible for a pion to produce a tachyon with the mass in the 
range required to fit the OPERA data. 

Another proposal to account for a superluminal particle has been suggested by Cole- 
man and Glashow (CG) S 13 • The idea is that the i-th particle has, in addition to its 
own mass m„ its own maximum attainable velocity c„ and obeys the standard dispersion 
relation: 

E^^pfcl + m^c^ . (5) 
The CG muon neutrino can indeed account for the OPERA data without any trouble 



associated with its production from pion, as we are going to discuss in section III To 
explain the observation of neutrinos associated in time with SN1987a, it is however 
necessary to introduce another neutrino with speed practically equal to c. This brings 
severe problems to neutrino oscillations, so that even the CG muon neutrino appears not 
to be a fully satisfactory explanation. 

We draw our conclusions is section [ 



II. PROBLEMS OF A TACHYONIC INTERPRETATION 



A. Energy independence of the early arrival times 

If the neutrinos produced at CERN are tachyons with mass m, after having travelled 
a distance L « 730 km, their associated early arrival time is 6t - ^ with ^ « 2.4 ms. 
Consider two tachyonic neutrino beams of energy £i and E2, with £1 < £2 for definiteness. 
As follows from eq.Q, the ratio of their early arrival times 6ti and 6t2 has a simple energy 
scaling: 




The early arrival time of a tachyon neutrino beam is indeed smaller the larger is its energy. 
In particular, for £2 ~ 3£i, one expects 6^2 ~ 6ti/9. The difference of the arrival times is 
thus negative: 6^2 - ~ -6^1. 

Now, the OPERA Collaborations considers two sample neutrino beams with mean 
energy equal to El = 13.9GeVand£2 = 42.9 GeV respectively ^. The ratio of these energies 



^ If neutrinos were tachyons, the energy reconstruction of the OPERA Collaboration should be revisited. 
However, for tachyonic masses smaller than GeV, such effect is negligible for the sake of the present 
considerations. 
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is indeed close to 3. However, the experimental values of the associated early arrival times 
are respectively dh = (53.1 + 18.8{stat) + 7A{sys)) ns and 6t2 = (67.1 + 18.2(stof) + 7A{sys)) 
ns. These data display no evidence of an energy dependence. OPERA quotes a value 
6t2 — 6h = (14.0 + 26.2) ns for the difference of the arrival times 6t2 — 6ti. Far from being 
close to -6^1 as expected for a tachyon, the latter value is even slightly positive, although 
consistent with zero. 

This simple argument disfavors the tachyon explanation of the OPERA data. The 
same conclusions were drawn in refs. HHHH (appeared when this paper was completed), 
carrying out a detailed numerical analysis and including in the fit not only the recent 
OPERA data but also the Fermilab data, which do not display any energy dependence 
too. 

One could however still question this conclusion, since the energies £i = 13.9 GeV 
and £2 = 42.9 GeV quoted by OPERA are mean ones and if we consider the 3a range 
associated to 6t2 - dti we find the interval [-65, 93] ns. 



Suppose that we close an eye on the energy dependence and we stick to the interpre- 
tation of the OPERA early arrival time in terms of a tachyonic muon neutrino. As we are 
going to discuss, arguments based solely on kinematics allow to obtain an indication for 
the value of the tachyonic muon neutrino mass. 

In terms of the muon neutrino energy E and the muon neutrino velocity v, the tachyonic 
muon neutrino mass m is simply given by (see eq.(|4])): 



Since OPERA deals with neutrinos with mean energy (E) = 17 GeV and observes (v-c) /c = 
(2.48 + 0.28 + 0.30) x 10"^, the corresponding tachyonic mass value is mc^ = (110 - 130) 
MeV at la, and (85 - 146) MeV at 3a (statistical and systematic errors are summed in 
quadrature)^ . This is also graphically shown in figjl] 

This range of values is close to the muon mass, m^,c^ « 105 MeV, and it would be 
fascinating to postulate that the muon neutrino mass is just the charged muon one, upon 
a rotation of n /2 in the complex plane. However, these conjectures are going to be severely 
challenged by other experimental data. 



^ Clearly this value is just a rough estimate, due to the significant energy spread of the neutrino beam, but 
we said that we ignore this as it would also cause an energy dependence of the early arrival time. 



B. Tachyon mass range from OPERA 




(7) 
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mc^ [MeV] 



FIG. 1: Relative speed of the tachyon with respect to light as a function of the tachyon mass and 
for selected values of the its energy: from left to right E = 13.9, 17.0,42.9 GeV. The shaded regions 
display the OPERA measurement, with 1,2,3a error bands (statistical and systematic errors are 
summed in quadrature). The vertical line marks the value of the muon mass, OTj,c^ ^ 105 MeV. 

C. Production from pion 

Indeed, first of all we must wonder whether a tachyonic muon neutrino with such a 
large mass could ever be produced in pion decay: n — > /iv^,. The OPERA muon neutrinos 
are in fact obtained from pions decaying in flight in a 1 km long vacuum tunneP. 

Let us focus on the kinematics of pion decay at rest. Clearly, momentum conservation 
requires the muon and the neutrino to be produced back to back and with the same 
momentum p. Energy conservation requires in addition: 

m„c^ = ^(pc)2 + (mf,c2)2 + ^(pc)2 - (mc2)2 , (8) 

where m^, ni^i stand for the pion and muon masses respectively. As in reaches its maximum 
allowed value when the muon neutrino has null energy (p = mc), one derives an upper 

bound in < ^jm^ - « 92 MeV/c^. This limit is marginally compatible with the 
tachyonic neutrino mass range derived before from OPERA results, see figjl] 

To better support this conclusion, in fig. |2] we display the /j and v^, energies and 
momenta as a function of the muon neutrino mass, assuming the latter to be a tachyon 



^ A proton beam of 400 GeV/c from SPS hits a graphite target producing pions which are focused by two 
magnetic horns and directed towards the tunnel. 
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FIG. 2: Kinematics for n ^.iv^i assuming that the muon neutrino is a tachyon (T) or a bradyon 
(B). The muon neutrino energy Ev^, (solid), the muon energy E^^ (dot-dashed) and the associated 
momentum (dashed) are shown as a function of the muon neutrino mass. 

(T) or a standard subluminal particle, also called bradyon (B). Deviations from the well 
tested values of the muon energy and momentum allow only a small tachyonic mass, say 
smaller than about 10 MeV/c^. When inserted in eq. (jij) and keeping the muon neutrino 
energy £ = 17 GeV, the bound m < 10 MeV/c^ would imply {v - c)/c < 1.7 x 10"'', which 
corresponds to an early arrival time at OPERA of 6t < 0.6 ns. This is at least two orders 
of magnitude below the observed value. 

We can rephrase all this also in another way: To end up with {v — c)/c « 2 x 10~^ 
while keeping m = 10 MeV/c^, a muon neutrino beam energy E « 1.4 GeV would have 
been necessary. The latter value seems to be definitely too small with respect to the 
reconstructed muon neutrino energy. We conclude that the tachyonic explanation of the 
early arrival times of muon neutrinos at OPERA is ruled out. 

We also note that a kinematical analysis of muon decay {[.i — > eVf.v^i) using the tachyonic 
mass range suggested by OPERA would produce serious difficulties. For simplicity, 
consider that v^, has a tachyonic mass m while e and Vg are massless. Then, in the corner of 
the phase space where £v,, = (and consequently pv,, = the electron energy can be as 
high as Eg < wj^,c^/2(1 + m/m^,) (intuitively, e and Vg have to balance the large momentum 
of the tachyon, and they have a lot of energy to do so, since the tachyon energy is zero). 
Even more significantly, one can see that for every value of the allowed energy range 
for the tachyon, < £v,, ^ wj^,c^/2(l - m^/m^J, the maximum value of Eg is larger than 
m^iC^I2{\ + np-jm^^). In conclusion, values of m > 10 MeV/c^ would be immediately 
detectable in the electron spectrum, whose endpoint would be much larger that m^/2. 
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Once again, this argument rules out the hypothesis that the OPERA muon neutrino is a 
tachyon. 



D. Supernova SN1987a requires an <eV electron-neutrino 



Let suppose that we close another eye on the problems associated with the production 
of a beam of tachyonic muon neutrinos with 100 MeV mass and persevere on this road. 
Can we agree with the SN1987a data? As we are going to discuss this is possible. 

The SN1987a is L = 1.68 x 10^ ly far from the Earth and exploded releasing a huge 
neutrino signal, with typical energies 10-20 MeV, which allowed the first direct detection 
of astrophysical neutrinos. All neutrino flavors were emitted but Kamiokande-II, 1MB 
and Baksan were designed to detect mainly electron anti-neutrinos. The signal lasted 
about 10 s and the photons also arrived within a few hours. See for instance ref. IITOl for 
a recent review and a list of references. 

The advance (delay) of a tachyonic (bradyonic) anti-neutrino with respect to light is 
6T - T\v - c\/c, where T - L/c is the time associated to the photon trip from SN1987a to 
the Earth^. The fact that photons and electron anti-neutrinos arrived without few hours 
implies that 6T/T - \v — c\/c 10~^, which in turn translates into an upper bound for the 
(bradyonic or tachyonic) mass of the electron anti-neutrino of about 1 keV. 

Electron anti-neutrinos arrived with a time spread AT < 10 s, as indicated by observa- 
tions. This poses a much tighter limit on their (tachyonic or bradyonic) mass m^^ than the 
one just discussed. The time spread AT = IT2 - Ti| of two neutrinos with energies £1 and 
E2 (with £1 < £2) is 

AT ml ( 

For the numerical values mentioned above, one obtains mp^, < 40 eV/c^. This limit, applies 
however to the electron anti-neutrino^. 

The SN emits all neutrino flavors. Let suppose that it emits also a 100 MeV tachyonic 
muon neutrino: its advance with respect to light would be of about 6t « Ayr, but with an 
enormous spread as can it can be realized by considering eq.(|9]) in the case of a particle 
with mass bigger than its energy. These neutrinos would have certainly escaped detection. 



* The SN exploded when the Earth was in the quaternary period and on such timescales effects due to the 

expansion of the universe can be neglected. 
^ We recall that an even stronger bound on the mass of a tachyonic electron neutrino follows from tritium 

beta decay experiments fill, which set a limit of few eV/c^. 
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E. Oscillations: game over 

According to the picture emerged so far, the ratio between the tachyonic mass of the 
muon neutrino suggested by OPERA and the mass of the electron anti-neutrino suggested 
by SN1987a would be as large as 10^. In principle, the formalism of neutrino oscillation in 
the tachyonic case is the same as for an ordinary neutrino [12], but it appears difficult to 
come to pact with the robust informations coming from neutrino oscillation experiments. 
Indeed, these experiments put stringent bounds on the difference of neutrino masses 
squared: lAmj^l « 2.4 x lO"'^ eV^ and Aml^ « 7.6 x 10"^ eV^ [13J. Even though the analysis 
of the experimental data should be redone since the fluxes at the production in the Sun 
for electron anti-neutrinos and in the atmosphere for muon neutrinos would change, it 
seems hopeless to find an agreement with experimental data on oscillations. 

III. MUON NEUTRINO A LA COLEMAN-GLASHOW 

As an alternative scenario, consider now two CG neutrino mass eigenstates with 
masses m-[ and m2 not larger than (9(eV)/c^, with different limit speeds ci and C2 UHIZl. 
One may infer |ci — c\/c < 10~^ from SN1987a as discussed in the previous section, while 
(C2 - c)/c « 2.5 X 10"^ as suggested by OPERA. 

Let us suppose that vi has a significant mixing with the muon neutrino of the OPERA 
beam, and that vi mixes significantly with the electron neutrino. We now revisit for this 
CG superluminal neutrino model the same issues discussed for the tachyon. 

First of all, the early arrival time of the muon neutrino beam is energy independent, 
since now C2 is a constant (already chosen to reproduce the results from OPERA) and we 
assumed m2C2 < 0{eV); with these assumptions the standard kinematics used for event 
reconstruction at OPERA need not be modified. 

At first glance, one could envisage no problem for the production of such CG muon 
neutrino from pion decay. However, a careful kinematical analysis reveals that the 
situation should be considered in more details. 

We first assume that c„ - c^i - c. For the sake of the comparison with OPERA, the 
relevant configuration is the one in which the CG muon neutrino and the muon have 
negligible transverse momenta with respect to the pion momentum, so that the CG muon 
neutrino is actually emitted in the Gran Sasso direction. In this case, there is an upper 
bound for the CG muon neutrino energy (see also [141 [15| ): 

(m^ - ml)c^ c 
£v< \^ ' ^«3GeV, (10) 

where the numerical value is obtained by considering ~ 60 GeV, together with the 
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OPERA result (C2 - c)/c « 2.5 x 10"^. Clearly, this bound is violated by OPERA, that 
detects muon neutrinos with energies much larger than this value. 

A widely different scenario follows if one assumes Cn — c and c^i, - ci- In this case 
the pion decay is forbidden unless the pion energy £„ is smaller than a threshold energy 
given by 



En < 



( (ml - 



1/2 



C2 - C 



13GeV, (11) 



where the OPERA result for (c2 - c)/c has been used. Clearly, muon neutrinos should 
have Ev ^ En- Again, this is in contrast with observation since the mean energy of the 
pions produced at CERN is about 60 GeV and OPERA detects neutrinos with energies up 
to about 80 GeV. 

This discussion shows that the observed phenomenology depends critically on the 
actual values of the c/s of the three particles involved in pion decay, whose values have 
been recently reviewed in |15|. Since at this stage a discussion of all possibilities would 
be inconclusive, we do not elaborate further on this point. 

As for the SN1987a, a beam of CG V2 would pose no problem, because it would have 
simply arrived about 4 yr in advance with respect to the photons and the other vi's. Most 
probably it would have escaped detection since the detectors had a lower sensitivity to 
muon neutrinos (moreover Kamiokande-II started taking data only in 1985). At variance 
with the tachyon case, it is important to remark that the V2 beam does not spread out in 
time but all the neutrinos arrive within a few seconds, because we assumed their mass to 
be smaller than eV/c^- 

A serious problem for CG neutrinos is again due to neutrino oscillations, as can be 
shown by using the formalism of refs.||6HZl- The two CG neutrino eigenstates travel at 
different speeds and this affects the neutrino oscillation probability similarly to a difference 
in mass: 



P(v,^v,) = l-sin^2gsin^^ r'"\:"^^^ + , (12) 




where is the mixing angle, R is the distance from source to detector, c - {ci + C2)/2, 
6c — C2 — Ci and E is the neutrino energy, typically in the range of a few MeV for reactor 
and solar experiments. For numerical estimates, it is perfectly safe to replace c with c. 
Oscillation experiments (see for instance [16J) indicate a value for rn^ — m^^ eV'^/c^. 
This translates in a sensitivity to 6c/c of about 10~^^, much smaller than what would 
be needed to explain the OPERA data. This can be seen as follows: the experimentally 



tolerated oscillation frequency is the one of the first term of the sin^ argument in eq.(12 
with [m^ — tn^c^ ~ 10~^ eV^ and E ~ MeV. A comparable frequency would result from the 
second term of the argument of sin^ only if bcjc ~ 10~^^. Also due to the different energy 
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dependence of these two terms, it seems unlikely that a cancellation might be at work for 
a much larger value of bcjc ~ 10~^ in the energy range probed by oscillation experiments. 
A similar analysis was done in refs.||6l|7l, suggesting an even tighter limit bcjc ~ 6 x 10~'^^. 

In conclusion, also CG superluminal neutrinos seem not to provide a fully satisfactory 
explanation of the OPERA results. 

IV. CONCLUSIONS 

The evidence for muonic neutrinos traveling slightly faster than light in vacuum as 
reported by the OPERA Collaboration, motivated us to explore two possible interpre- 
tations of the data: the hypothesis that the muon neutrino is a tachyon or that it is a 
Coleman-Glashow neutrino. 

We demonstrated that the tachyonic interpretation is hardly reconciled with the energy 
independence of the OPERA data, as shown also by (H |9l. The real problem that we 
point out here is that it would be impossible to produce a 100 MeV tachyon from pion 
decay. The data associated with SN1987a can be interpreted by assuming an eV electron 
anti-neutrino. This picture however clashes with what is known concerning neutrino 
oscillations. 

A Coleman-Glashow superluminal neutrino beam would presumably face problems 
with kinematics, but it is difficult to asses details here because the kinematic bounds 
strongly depend on the actual values of the limiting speeds of the particles involved in 
pion decay. A CG neutrino could be easily reconciled with SN1987a data but, on the other 
hand, it would be not possible to reconcile the model with neutrino oscillations. 

In conclusion, the picture emerging from combining OPERA with other experimental 
data is that neutrinos should not obey a tachyon type nor a Coleman-Glashow type dis- 
persion relation, but rather a dispersion relation with a very peculiar energy dependence 
Il9l[l7l. Even in this case, a serious problem could be represented by energy losses due to 
electron positron pair production [,18||. 
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